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ON THE POISSON EQUATION AND DIFFUSION 
APPROXIMATION 3 

By E. Pardoux and A. Yu. Veretennikov 

Universite de Provence and University of Leeds 

We study the Poisson equation Lu + / = in R'', where L is the 
infinitesimal generator of a diffusion process. In this paper, we allow 
the second-order part of the generator L to be degenerate, provided a 
local condition of Doeblin type is satisfied, so that, if we also assume a 
condition on the drift which implies recurrence, the diffusion process 
is ergodic. The equation is understood in a weak sense. Our results 
are then applied to diffusion approximation. 

1. Introduction. This is the third in a series of papers devoted to the 
study of the Poisson equation in and diffusion approximation. In this 
paper we consider the degenerate case. 

The study of diffusion approximation [i.e., obtaining the hmit of 
in (19)] was initiated by Khasminskii [5], and developed by many authors, 
including Papanicolaou, Stroock and Varadhan [10] and Kushner [8]. Such 
results, and the formulation of the limiting stochastic differential equation, 
require the solution of a Poisson equation Lu + / = 0, where L is the in- 
finitesimal generator of a Markov process (at least in the case where the dis- 
turbance is Markovian; in the non-Markov case a substitute of the Poisson 
equation replaces it), whose right-hand side / is the highly oscillating coeffi- 
cient of the approximating differential system. When the disturbance in the 
approximating ODE is compact valued, the Poisson equation is formulated 
in a compact set, and the corresponding theory is well known; the result can 
be proved under quite explicit conditions on the coefficients (see [2], Chap- 
ter 12, Section 2). When, however, the disturbance of the diffusion takes 
values in all of W^, there was until recently no way of deriving estimates for 
the solution of the Poisson equation in terms of explicit conditions on the 
data; see Chapter 12, Section 3 in [2]. 
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This was the starting point for our work. We focused on the case where 
the disturbance is an ergodic M'^-valued diffusion process, and the ergodicity 
fohows from expUcit conditions on the coefficients. In the first paper [13], 
using mainly probabihstic arguments (together with some estimates from 
the theory of partial differential equations) , we solved the Poisson equation 
in R'^ for the generator of an elliptic and ergodic diffusion, and obtained 
estimates (which we believe are rather sharp) of the solution. We then used 
that result to establish a diffusion approximation result under very explicit 
conditions on the coefficients. In the second paper [14], we considered the 
case where the coefficients of the equation (18) depend on the process. 
This forced us to study a Poisson equation where both the PDE operator 
and the right-hand side depend on a parameter, and establish regularity 
results of the solution in terms of that parameter. We were forced for that 
purpose to use essentially results from the PDE theory. 

The aim of this third paper is to consider the situation of the first pa- 
per, where we now relax the ellipticity assumption. While a condition on 
the behavior of the drift at infinity [condition (Ai,)] implies the positive re- 
currence, irreducibility, which was in our previous works a consequence of 
ellipticity, is now a consequence from a type of "local Doeblin condition" 
[condition (Di)]. While those conditions are not explicit conditions on the 
coefficients of the diffusion, they are implied both by the ellipticity assump- 
tion and by the "restricted Hormander condition" (i.e., the assumptions that 
the diffusion vector fields, together with their brackets of arbitrary order, 
span the whole space at each point). We further give one example where none 
of these conditions holds, while our condition (Di) is satisfied, together with 
the additional "regularity" condition (At)- 

We then consider a weak formulation of the Poisson equation, which is 
solved by the same probabilistic formula as in the elliptic case. We finally 
apply those results to the diffusion approximation problem. We prove weak 
convergence in the sense of the 5-topology of Jakubowski [4]. The difficulty 
in proving convergence in a stronger sense is related to the lack of smooth- 
ness of the solution of the Poisson equation. Let us also mention that our 
diffusion approximation is "less general" than the one considered in [13] 
(except that we relax the ellipticity condition, as explained above), in that 
the approximating differential system does not contain a stochastic integral 
(i.e., the coefficient H in [13] does not appear here). 

Let us point out that precise regularity of the solution of the Poisson 
equation under the Hormander condition follows from Theorem 18 in [16]. 

The paper is organized as follows. Section 2 contains our assumptions, as 
well as some essential ergodicity results from Veretennikov [23]. Section 3 is 
devoted to the study of one example of degenerate coefficients, which satisfy 
our assumptions. The Poisson equation is studied in Section 4, while the 
diffusion approximation result is derived in Section 5. 
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2. Moment bounds and convergence to the invariant measure. Consider 
the stochastic Ito equation 

(1) dXt = b{Xt)dt + aiXt)dBt, Xo = x£R'^, 

where {Bt,t > 0} is a A;-dimensional Brownian motion, 6 is a locally Lips- 
chitz vector-function of dimension d and a is a d x k matrix-valued locally 
Lipschitz function. We assume that aa* is bounded and possibly degenerate, 
and that the unique solution of (1) satisfies 

(At) \/R>0 sup E^.inf{t>0:|Xt|>i?+l}<cx). 

\x\<R 

Let us introduce the following recurrence condition: 
(Af,) lim {b{x),x) = —oo. 

Note that this condition prevents the solution of the SDE (1) from explod- 
ing, so that the process {Xt} is well defined for all t > 0. Let R> and 
T = r^ = inf(t >0:|Xt| <R). 

Finally, we assume the following "local Doeblin" type condition. 

Let BcW^ and 

:=M{t>0:Xt£B}, 

and (in the following formula, ts > depends only on B) 

ri^i := inf{t > rf + : G B}. 

Define the ^^process in i?" in discrete time as X^ := X^b . Denote by P^{n, x, dx') 
the n-step transition probability of {X^). We say that the local Doeblin con- 
dition holds true for the process {Xt} if for any R' > there exists R> R' 
such that the process in i? = B/j := {x G M"' : \x\ < R} satisfies the following: 
there exists an integer uq = no(R) > such that 

(De) inf / mml^^P^^^^i^,l\p''{no,x',dx'')=:q{R,no)>0, 
^ ' \x\,\x'\<rJb \P^{no,x',dx")' J ^ ' ^ '^v . u; 

where ^ ("■O'^''^^ ) jg defined as follows. Let 

F" {no,x' ,dx") 

P^{no,x,dx") = ipx^x'{x")P^{no,x',dx") + i^x,x'idx") 

be the decomposition of P^ {n(),x,dx") into its absolutely continuous part 
w.r.t. 

pB 

(no, x', dx"), and the part Vnc,x'{dx") which is singular w.r.t. P^{no, x', dx"). 

Then 

P^{no,x,dx") _ „ 
P^{no,x\dx") -'^-'-'^'^ 
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The assumption (Di) requires, in particular, that the mass of the singular 
part is not close to 1, and moreover, it imposes a certain quantitative esti- 
mate on the total variation norm for the difference of two measures uniformly 
on the compact B. We assume throughout the paper that the process {Xt} 
satisfies this assumption. We shall give in the next section one example with 
a nonelliptic diffusion coefficient, for which conditions {At) and {Di) hold. 
The proof that example satisfies our conditions will use stronger conditions, 
which are easier to verify. 

We note that the two assumptions {Af,) and (Df) imply the existence and 
uniqueness of an invariant probability measure. For this and the proof of the 
next proposition, see [23]. Note that in [23] the "process in B" is defined in a 
slightly different manner, since it is extracted from the sequence {Xn,n £ N} 
rather than from {Xt,t > 0} as defined here. However, the adaptation of 
those proofs is rather obvious. 

Proposition 1. Under the assumptions (At), (At,) and (Di), for all 
m' > m + 2 > 2, there exists C such that for all x G W^, t>0, 

(2) E^\xtr <c{i+\xr'). 

Moreover, 

(3) K^\Xtr <oo Vm>0, 
and for any k > 0, 2k + 2 <m, 

(4) var(/if - ^) < C(l + |x|™)(l + 

where "var" denotes the total variation norm of a signed measure over the 
Borel a-field, fif is the law of Xt when Xq = x, fi is the unique invariant 
measure of X andE,^ means the expectation w.r.t. fi. 

Proposition 2. Let the assumptions (At), (Ah) and (Di) be satisfied. 
Then for any p> 0, 

'Exi sup \Xti\^ \ = o{Vi) ast^oo. 
\o<t'<t J 

The proof of Proposition 2 is similar to that in [13]; hence, we drop it. 
The following corollary will be used in Section 4, for the proof of tightness. 

Corollary 1. Under the same assumptions, for any T > 0, p> 0, 

eKx i sup \Xt/i;2 [^ ) ^ as £ ^0. 

\0<t<T J 
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3. Sufficient conditions and one example. In this section, we first state 
two conditions, wliich we prove to be, respectively, stronger than (At), and 
stronger than (Di). Then we give one example with a degenerate diffusion 
coefficient, which satisfies those stronger conditions. 

Proposition 3. A sufficient condition for condition (At ) to hold is that 
for each R>0, there exists f G C{R'^, [0, 1]) with supp(/) C {x; \x\>R + 1} 
and to such that 

inf EJ{Xt,)>0. 
\x\<R 

Proof. Let i? > 0. It follows from our condition that there exists c, to > 
suchthat for allx eM'', \x\ < R, Prd-'^tol > ^+1) > c; this implies > 
R\) >c, and the same is true for R+1 instead of R, with new c and to- Let 

SR+i = mi{t;\Xt\>R + l}. 

It follows from the Markov property of {Xt, t>0} and the previous estimate 
that for all a; G M'^, 

P,(5ij+i>nto)<(l-cr. 

The result follows. □ 



We now formulate what we call the condition (Dgi) ("strong local Doeblin 
condition"). 

For each i? > 0, there exists Ar C Br, tn > 0, c{R) > 1 such that for 
all X G Br, the transition probability of our diffusion process {Xt,t > 0} 
satisfies 

p{tR, x; dy) = q{tR, x, y)^{dy) + u{tR, x; dy), 
(Dsi) ^ 

-7^<q{t,x,y)<c{R), y € Ar, 
c{R) 

where // is a probability measure on M*^ such that //(A/j) > 0. 

Note that the upper bound here is not actually necessary for our aims; 
however, it holds along the lower bound in all cases known to the authors. 

Proposition 4. The strong local Doeblin condition (Dg^) implies the 
local Doeblin condition (Di). 

Proof. We choose an arbitrary > 0, and denote B = Br. We decom- 
pose the transition probability of the process in B (defined with t_B = tR) as 
follows. For X £ B, 

P^{l,x,dx')=¥^{Xt^ G dx',Xt^ G B) + i^'{x,dx'). 
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It follows from (Dgi) that 

c[K) 

Hence (Di) holds with no = 1 and q{R, 1) = □ 

Example 1. Let b and cq satisfy the above conditions (Ax) and (Ah), 
b e C^(M'^,M'^), fJo G C2(]R'^,R'^x'^), and let ctq be uniformly nondegenerate. 
Let a : M'^ — > [0, 1] be a mapping, such that the set {a = 0} is the union 
of countably many disjoint connected closed subsets of Mf^, such that each 
bounded subset of M.'^ intersects at most finitely many of those, and the 
set {a > 0} is connected. We now assume that for some 5 > such that 
{a > 6} is connected and for each R' > there exists R> R' such that the 
set {|x| = R} does not intersect the set {a < 5}, and moreover that for any 
i? > 0, there exists M such that the solution of 

^{t) = b{x{t)) 

exits in time less than M from {a < 5}, whenever x{0) £ {a < 5} D Br. Let 
(j{x) = a{x)aQ{x). Then the pair {b, a) satisfies the assumptions (At) and {Di 

We first prove the following. 

Lemma 1. The condition of Proposition 3 is satisfied in Example 1. 

Proof. We consider the stochastic equation for the process {Xt,t > 0}, 
written in Stratonovich form (the reason for this is that we shall soon use 
Stroock and Varadhan's support theorem), that is, 

dXt = b{Xt)dt + aiXt)odBt, 

where 

bi{x) = bi{x) - ^ 

We will use below the notation (V(7)(7o for the vector 2a~^{b — b). We now 
consider the controlled ODE 

^(t)=%(t)) + a(y(t)Mt), 

y(0)=Xo, 

where we choose the feedback control u{t) = $(y(t)), with 

\aQ^{Va)aQ{x), if ^(x) > 0, 



<I)(x) 



if a{x) = 0. 
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It is easy to check that {y{t), t > 0} coincides with the solution of the ODE 
-{t) = b{x{t)), x{0)=X{0). 

Let 

r = t(x(0)) = inf{t > 0; a{x{t)) > 5}. 

Choose p = (2|| Va||oo,fl)"^5, g G C(M'^; [0, 1]) with supp(5r) C {x, \x\ < p} and 
5( = 1 on a neighborhood of 0. We have used the notation 

||Va||oo,fl = sup ||Va(x)||. 

\x\<R 

The above considerations and Stroock and Varadhan's support theorem 
(cf. [17]) imply that 

lEx(0)5(^r-x(T))>0. 

Moreover, that last quantity depends continuously on x(0), hence it is bounded 
away from zero for x{0) £ {a <5} D Br. Hence our construction yields that 
with a probability which is bounded away from zero, a{XT-) > S/2, where 
r < M is a deterministic time which depends only on 2;(0). r = whenever 
a(x(0)) >6. 

Let / e C(M'^, [0, 1]) satisfy supp(/) C {x; \x\ > R + 1} and f{x) = 1, when- 
ever |x| > i? + 2. Using again Stroock and Varadhan's support theorem, we 
have that 

inf inf ^ , Ex'/(^2M-t(x)) > 0. 

Proposition 3 with to = 2M and the above / now follows from the Markov 
property. □ 

Lemma 2. The pair {b,a) from Example 1 satisfies the condition (Dgi). 

Proof. It follows from the proof of Lemma 1 that there exists ^ > 
and a mapping r G C{Br, [0,M]) such that for all x G Br, 

P,(q(X,(,))> V2)>e. 

Next, we choose a closed ball A C vat Br n {a > 0}. Using again Stroock 
and Varadhan's support theorem, we deduce that there exists N > M such 
that 

inf inf F^,{Xj^_^(^)£A)>0. 

Combining the above two statements with the help of the Markov property, 
we obtain that 

(5) inf P^(X7vG^)>0. 
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Next we choose another closed ball A' such that A C int A' C A' C Br n 
{a>0}. 

For any function 93 G C(A',M+), with supp(v3) C A, we consider the solu- 
tion {tt(t,x),0 <t< l,x £ A'} of the backward linear parabolic PDE [here 
a{x) = aa*{x)] 

- (t , X) + - 5: a., (x) (M) + 5: 6, (.) — (t , X) = 0, 

< t < l,x G A', 

u{l,x) = ip{x), u{t,x) = 0, x£dA'. 

We have that for all x G A' , n(0,x) = Exip{Yi), where the process {Yt,0 < 
t < 1} is the solution of the SDE (1), which is killed when it reaches the 
boundary of the set A' . It follows from the parabolic Harnack inequality 
(see, e.g., [7], page 131) that there exists > such that 

x,x'eA u[0,x') 

that is, 

sup 7T7T<iV, 

x,x'eA ^x'^iyi) 

for all (f G C(A',M+), with supp((/7) C A. We choose one particular point 
xq G a, and define fj.{dy) = P^oC^i ^ dy). It follows from the above that for 
each X G A, Fx{Yi G dy) is absolutely continuous with respect to /U, and 
moreover the Radon-Nikodym derivative q{x,y) satisfies 

(6) N~^<qix,y)<N, 

for all x,y £ A. 

Condition (Dgi) now follows from (5), (6) and the Markov property. □ 

Remark 1. It is rather clear that one can verify our assumptions in 
many other situations, where det(a(x)) may vanish in a similar fashion as 
a(x) does in Example 1. All that is to be verified is a condition like (5), 
both for a set of the same type as A and for -Bjj+i- 

Remark 2. In the strictly elliptic case the same arguments based on 
Harnack's inequality establish the condition (Dgi), provided a = aa* /2 is 
continuous, and b locally bounded. The same is true, with fj, = Lebesgue 
measure, whenever the coefficients are smooth, and the Lie algebra of vectors 
fields generated by the columns of the matrix a has full rank at any point 
of R'^. 



POISSON EQUATION AND DIFFUSION APPROXIMATION 9 

4. The Poisson equation in M'^. We consider the Poisson equation in M*^ 

(7) Luix) = -fix), 
where 

L = Y^ a-ij ix)dx, + hi {x)dxi ^ 

with 

a{x) = aa*{x)/2, 

and / e C{W^) satisfies 

<C(l + |x|)^ for some /3 e M, 

so that due to Proposition 1, / is integrable with respect to the invariant 
measure fi, and 

(Af) J f{x)fi{dx) = 0. 

In the nondegenerate case, the solution of (7) has the representation 

(8) u{x)= EJiXs) ds. 

Jo 

In the degenerate case it is useful to extend the notion of equation (or 
solution; we prefer the former): we say that u solves the integral Poisson 
equation if for any t> 0, x G W^, u{Xt) is P^-integrable and 

(9) u{x) = KM^t) + t Exf{Xs) ds. 

Jo 

This notion is similar to probabilistic or martingale solution of a parabolic 
equation in [18], also for the degenerate case; in this respect it is worth re- 
membering that a classical solution to the degenerate parabolic equation 
was first constructed by Gikhman [3] . It is also easy to show that a continu- 
ous function solution of the integral Poisson equation is a viscosity solution 
of the Poisson equation, in the sense of [1]. 

Notice that (9) may be reformulated in the following form: for all x G M'^, 

(10) u{Xt) — u{x) + / f{Xs) ds is a martingale under P^;. 

JO 

Indeed, (10) implies (9) by taking expectation. Vice versa, if we substitute 
zero by t' , and x by Xti in (9) {t' < t), then by virtue of the Markov property 
we get 

u{Xt') = Ex^MXt) + Ex,, £ f{Xs) ds, 
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or 



E 



-u{Xt')+u{Xt)+ f'f{X,)ds 
Jt' 



Ft' 



0. 



Hence, it follows that 
E 



u{x) + u{Xt)+ ( f{Xs)dsFt' =u{Xt,)-u{x)+ f f{X,)ds, 
Jo J Jo 



which means exactly the desired martingale property. 
Define 

POO 

u{x)= / \E,f{Xt)\dt. 
Jo 



Theorem 1. Let the assumptions (At), (Af,) and (D^) he satisfied. We 
assume that there exists < /3 such that \f{x)\ < C{1 + \x\^) with C > 1 and 
that (Af) holds true. Then (8) defines a continuous function u, which is a 
solution of (9) and satisfies the following properties. For any m > /3 + 4, 
there exists Cm which depends only on m,(3, the value supj ^. and on 

the constants C in (2), such that 

(11) \u{x)\<u{x)<Crr,{l + \xD, X G M'^, 

SO that in particular u is fi-integrable. Moreover, again for any m> (3 + A, 



(12) sup(l + \x 



\m\ — l 



u{x) 



N 



Exf{Xt)dt 







as N ^ oo. 



In addition, u is centered in the sense that 



(13) 



u{x)fi[dx) = 0. 



The solution is unique in the class of solutions of (9) which satisfy properties 
(11) and (13). 

Theorem 2. Let the assumptions of Theorem 1 he in force. 
(i) // there exists C such that 

(14) |/(x)|<C(l + |x|)^-2 
for some /3 < 0, then u is hounded. Moreover, 

(15) sup|n(x)| < C7sup[|/(x)|(l + \x\)-'^+% 

X X 

where the constant C depends only on the constants C, m, k from (2)-(4) 
in Proposition 1. 
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(ii) // there exist C, (3 > with 
(16) \f{x)\<C{l + \x\f-\ 

then there exists C such that 

\u{x)\<C'{l + \x\)f^. 

Moreover, 



(17) 



1 + 1x1/3-2' 



where the constant C" depends only on the constants C , m, k from (2)-(4) 
in Proposition 1. 

The assertion of Theorem 1 is used in Theorem 2, which means that 
the last theorem gives additional information under additional assumptions. 
Theorem 2 gives in particular a criterion for u to be bounded. 



Proof of Theorem 1. The calculations are similar to those in [14]; 
however, they are not identical. Therefore we present the proof for the 
reader's convenience. 

A. u is well defined and satisfies (11). This follows from [22]; see Propo- 
sition 1. Indeed, 



u{x)= \E^,f{Xt)\dt 
Jo 



dt 



fiyMidy)-Kdy)] 



dt. 



Without loss of generality, we assume that (3 + 2 <m. Due to the inequalities 
in Proposition 1, one can choose p> 1, q> 1 with + = 1, such that 
pP < m and q<k+l. 

Indeed, if /3 = 0, then it is evident. Consider the case /? > 0. Let p = m/ (3. 
Then = I — f3/m, and {k-'rl)/q > 1 is equivalent to (/c + 1)(1 — /3/m) > 1. 
Since k + 1 is an arbitrary number less than m/2, then the last inequality 
can be satisfied if (m/2)(l — [3/m) > 1, which is equivalent to m > /3 + 2, 
and this is our assumption. Now, using Holder's inequality, and denoting all 
new constants by C (they may be different on each line), one has 



f{yWt{dy)-i^{dy)] 



dt 
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1/p/ r \V9 



\f{y)n^l-^{dy)+^^{dy)]^ (^j - ^l\{dy)^ dt 
<C + \yrM{dy) + /i(dy)]) '^"(varC^f - dt 

<c / (i+E^.|Xtr+E.|Xtr)i/p((i+|xr)(i+t)-('=+^))^/^dt 

Jo 

< c(i + ixry/'' /~(E^|Xtr + 1)^/^(1 + ty'-'^+'^^/ut 

Jo 

<c{i + \xr'r'' / {E^\xtr)^/p{i+t)-^^+^'^/ut+c{i+\xr') 

Jo 



1/g 



<C(l + |xr )• 

Thus, u is locally bounded and, moreover, (11) holds true with any m' > 
/3 + 4. The assertion (12) follows from the same calculations with instead 

B. u satisfies (13). Notice that if some function g is integrable w.r.t. the 
invariant measure fi, then for any s > 



E^[g{Xs)]n{dx) = J g{x)fi{dx). 
Due to (11), the function u is ;U-integrable. So, by virtue of Fubini's theorem, 
j E,f{Xs)dsfi{dx)= J E,f{Xs)fi{dx)ds. 

But clearly 

' EJ{Xs)fi{dx)= I f{x)fi{dx) = 0. 



C. u is continuous. It follows from the locally uniform convergence (12). 

D. u solves the integral Poisson equation (9). Let t > be a nonrandom 
value. First note that 

u{x)= / E,/(X,)ds+ / E^J{Xs)ds, 
Jo Jt 

where both integrals are well defined. On the other hand, from the Markov 
property of X, 

roo rOG 

/ EJ{Xs)ds= / E^Ex J{X,)ds 
Jt Jo 

pN 

= lim / E^ExJ{Xs)ds 

rN 

= lim E, / ExJ{Xs)ds 
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= lim E,n^(Xi) 
= E,u(^t), 
where u^{x) := E,xf{Xt) dt. Hence, 

u{x)-EMXt)= f'E^f{Xs)ds. 
Jo 

This is exactly (9). 

E. Uniqueness. For the difference of two solutions, v = u — u' , we have 
due to (9), v{x)=K^v{Xt). So 

v{x) ='Exv{Xt) ^ / v{x)iJ,{dx) = 0, t^oo. 

Hence, v{x) = 0. □ 

Proof of Theorem 2. The proof is identical to that in [14]; in partic- 
ular, the strong Markov property of the process Xt makes possible the use 
of the formula 



n(x)=E^n(X^fl)+E^ r f{Xt)dt, 

Jo 



which leads to boundedness condition for the function u in the first assertion. 
We refer to the calculations in [14]. □ 

5. Diffusion approximation. Let {Xt,t > 0} denote the solution of the 
SDE 

dXt = b{Xt)dt + a{Xt)dBt, Xo = x, 

and define Xf := Xi/^2, t > 0. Note that for some Brownian motion {B^} 
depending on e, solves the SDE 

(18) dXf = e~\Xf) dt + e~ V(Xf ) dBf, X^ = x. 

In this section, we are going to apply Theorem 1 to the singularly perturbed 
ODE 

(19) ^ = F{Xf , y/) + e-^Xf ,Yf), 0<t<T, Y,^ = y. 

Here e is a small parameter. The process X is the same as that of the 
previous sections, and we will again assume the same conditions (At,), (At) 
and (-D^). F and G are Borel vector-functions. The dimension of X is again 
d, the dimension of Y is i. We denote again by L the generator of the 
process X . The problem we are interested in is the weak convergence of 
the slow component Y^ as e— > 0. Concerning (19), we require the Lipschitz 
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condition with respect to the variable y, with a constant which may depend 
on x: 

(Al) \F{x,y)-F{x,y')\ + \G{x,y)-Gix,y')\<C{x)\y-y'\, 

where x — > G{x) is locally bounded. We now assume that for all x G M'^, 
G{x, •) G C^(]R^;M^), that dyG G C(M'^+^;M^') and the functions F, G satisfy 
the following polynomial growth conditions: 

\F{x,y)\<K{l + \y\){l + \xn, 

{Ap) \G{x,y)\<Kil + \y\)il + \xn, 

\\VyG{x,y)\\<K{l + \x\'^^). 

We assume moreover that for all y G and j = 1,2, . . . 

(Ac) [Gj{x,y)fiidx) = 0, 



where fx{dx) again denotes the (unique) invariant measure of X. It then 
follows from Theorem 1 that the Poisson equations 



LGj{x,y) = -Gj{x,y), j = !,■■■ 



which we in fact interpret as integral Poisson equations, see (9), have unique 
centered solutions 

G,{x,y)= E,Gj{Xl,y)dt. 
Jo 

Moreover, for some K and q'^, the following holds: 

\G{x,y)\<K{l + \y\){l + \x\''2), 



(20) 



VyGix,y)\\<Kil + \x\ 



93 > 



The values of q2 and ^3 can be deduced from those of (72 and qs by using 
Theorem 1 or Theorem 2, and the fact that Vj^G = VyG. 

In the next theorem, we make use of the S topology of Jakubowski on the 
space D([0,r];M^) of "cadlag" M^-valued functions defined on [0,r]. We refer 
to [4] for a definition of that topology and the presentation of its properties. 



Theorem 3. Let the assumptions (Ah), (Di), (At) and (Al), (Ap), 
(Ag ) he satisfied. Then for any T > 0, the family of processes , < t < 
r}o<e<i is uniformly S -tight in D([0,r];M^). IfY is an accumulation point 
of the family {V^ ,e 0}, then it is a.s. continuous, and it is a solution of 
the martingale problem associated to the operator 
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where 

b{y) = F{y)+Y^ J Gi{x,y)dyfi{x,y)i,{dx) 

i 

with 

F{y) = J F{x,y)iJi{dx) 

and 

d{y) = J [Gix,y)G*ix,y) + G{x,y)G*{x,y)]n{dx). 

If, moreover, the martingale problem associated to C is well posed (it is easy 
to state sufficient conditions for that), then ^ Y in the sense of the 
S -topology, and Y is the unique (in law) diffusion process with generator C. 

Notice that all integrals in the definition of C are well defined, as follows 
from Proposition 1. 

Proof of Theorem 3. 

Step 1. Preliminary computation. Let / G Cp(M^) (the set of functions 
of class which, together with their partial derivatives of order 1, 2 and 3, 
have at most polynomial growth of some order) and define 

f{x,y) = f{y) +eu{x,y), 

where eu{x,y) is a corrector to /, defined as follows, u is the solution of the 
Poisson equation 

Lu{x,y) = -{Vyf{y),G{x,y)), 

or in other words 

(21) u{x,y) = {Vyf{y),G{x,y)), 

where G : M"^ x ^ solves 

LG{x,y) = -G{x,y) 
in the integral form (9). Note that 

J dyG{x,y)fi{dx)=0, yeR^, 

and 

dyG{x,y) = dyG{x,y). 
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For each 6 > 0, we associate a mesh = to < < ■ " ■ < ^n. < • • • ) such that 
U — ij-i < (5, i > 0, and ti — > oo, as i ^ oo. For each t > 0, let N{t) denote 
that smallest integer such that t < tAr(t). It follows from our definition of the 
Poisson equation solved by G{x,y) that for all e > 0, (5 > 0, the following is 
a local martingale: 



e,5 



i<N{t)~l 



+ 



{Vf{m,G{Xl,Yn)ds 



Moreover, 



i<N{t)~l 

i<N{t)-2 

and for i < N{t) - 2, 

eu{Xl_^^,Y,\^^) - eu{Xl_^^,Y,\) 



Finally, 



e r^\Vyu{Xl^^,Y!),F{XlY!))ds 

J ti 

J ti 



■-eu{Xt,n^^^^J-Eu{X^,,Y^) 
-e E t^\^yu{Xl^,X),F{Xl.Y!))ds 

i<Af(t)-2 



{Vyu{Xl^^,Y!),G{X!,Y!))ds 



1 V- 



ti+iAt 



{Vyf{YfJ,G{Xl,YQ)ds. 



We now let 5^0 in the last identity, from which we deduce that the fol- 
lowing is a local martingale: 

Mf = euiXlYf)-euiX^o,yo) 

-e l\Vyu{XlY!),F{XlY!))ds 
Jo 
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- [\vyu{X!,Yn,G{Xl,Y!))ds 
Jo 

+ 7 [\vyf{Y!),G{XlY,^))ds. 
e JO 

Moreover, we have that 

/(y/) = f{Yon + f^(vf{Y!lF{XlY!) + -G{XIY!)^ ds, 

hence 

JO 

+ f\d'f{Y,')G{Xl,Y!),G{Xl,Y,'))ds 

JO 

+ Mf '-^ 

(22) 

+ e(V,/(yo^), G(Xo^ yo^)) - e(V,/(y/), , y/)) 

+ e A(v/(y/),v,G(xf,y/)F(xf,y/)) 

JO 

+ (92/(>;^)G(Xf,y/),F(Xf,y/))]cis, 

where {Mf'-^,t > 0} is a continuous local martingale which is localized by 
the sequence of stopping times 

5^:=inf{t;|y/|>n}, n = l,2,.... 

Step 2. 5-tightness. We shall make use of the iS-topology on ©([O, T]; M^), 
introduced by Jakubowski. The following result is a consequence of the re- 
sults in [4] and [9]: 

Proposition 5. The collection {Yf,0 <t< T}{o<£<i} is uniformly S- 
tight if it satisfies the two conditions: 

(i) For all 5 > 0, there exists M > s.t. 

p( sup \Y^^\>m]<5, 0<e<l. 
\o<t<r / 

(ii) Yf - Y^ = Ef + Vf + Aff , with 

(23) Ef ^ in probability, uniformly for t £ [0,T], 
and for each n G N, 

(24) sup E(||y^||TA5ji + {M')tasO < 

0<£<1 
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where WV^Wt denotes the total variation of between and t, and (Af^^ 
denotes the quadratic variation of the continuous local martingale 



We first prove that the sequence (X-^) satisfies (i). For that sake, we will 
use (22), with the function f{y) = log(l + |yp). Recall that the function u 
depends on /. Notice that for this choice of / one has 

(1 + \y\)\dyf{y)\ + (1 + \y\f\\dlf{y)\\ + (l + \y\f\\dlf{y)\\ < C, 

and then in particular [see (20)] 

\u{x,y)\ < K{1 + \x\'^3). 

Consequently, the absolute values of the integrands in the right-hand side 
of (22) do not exceed C(l + I'') with some q < oo. So {m/'^} is in fact a 
martingale, and there exist two constants C and q such that for < e < 1, 



E 



sup iog(i + |y/ 

0<t<T 



e|2^ 



<C sup E(l + |Xf|'?) <oo. 

0<t<T 



This implies that the condition (i) in Proposition 5 is satisfied. 

It remains to prove that (ii) is satisfied. For that sake, we choose f{y) =y 
in (22), yielding 

(25) Yf = Yo' + E! + Vf + M^, 

where 

E! = eG{X^,Y,')-eG{Xt,Yf), 

Vf= l\l + eVyG{Xl,Y!))F{XlY!)ds+ [\yG{XlY!)G{Xl,Y!)ds, 
Jo Jo 

and {Mf,t > 0} is a continuous local martingale. 

Now (23) follows from Corollary 1, (20) and (i), and the first half of (24) 
follows from {Ap) and (20), and we finally compute (M^). 

From (22) with /(y) = |yp, 

\Yf\' = \Y,'\' + 2 l\Y!,F{XlX)+'^yG{XlX)G{XlY!))ds 
Jo 

+ 2 j\G{XlY!),G{Xl,Y!))ds 
Jo 

+ Mt'^ 

+ 2e{Y,^G{X^„Y,n) - 2e{YfG{X!,Yf)) 

+ 26 A(y/,v,G'(xf,y/)G(xf,y/)) + (G(xf,y/),G(xf,y/))]ds, 

Jo 
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where {Mt',t>0} is a continuous local martingale. 

Now from Ito's formula for continuous semimartingales and (25), we de- 
duce that 

\Yf + eGiXlYf)\' 

+ 2 /V/, F{X!,Y,^) + VyG{XlY!)G{XlY[)) ds 
Jo 

+ 2e l\Y!,VyG{Xl,Y!)F{Xl,Y!))ds 
Jo 

+ 2e [\g{XI,Y,'), 
Jo 

{I + eV.GlXf , y/))F(Xf , y/) + VyG{X!,Y[)G{XlY:)) ds 
+ 2 / y/dMf + (M")j, 

JO 

where y/ = y/ + eG'(Xf ,y/). Comparing the last two identities, we deduce 
that 

(M-)i = 2 r(G(xf ,y/),G(Xf ,y/)) ds 
Jo 

+ e'\GiXf,Yf)\'-e'\G{X^„Y,^)\' 

+ e j\e {XI , y/) ds + Mf - 2 y/ dMi , 
Jo Jo 

where 

|V'.(rE,y)|<C(l + e)(l + |y|2)(l + |x|35). 

The second half of (24) now follows from (20) and the assumptions on the 
growth of G. 

Step 3. Identification of the limit. Let < s < t <T, and let be a 
bounded and S-continuous functional defined on D([0,r];M^), which is mea- 
surable with respect to the c-algebra a{x{r),x E D([0,T];R^) < r < s). Let 
/ G C^(]R^) be a smooth function with compact support. It follows from (22) 
that for all o > 0, such that t + a<T, 

E{[f{Y,%J-f{Y,%,)mY)) 

= EfcD,(y) r\vf{Y,'),F{X',,Y,') + VyGiX^.X)G{X'„Y,'))dr) 

\ Js+a J 
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+ E('^siY) r\d^fiY,')GiX^.,Yn,G{Xf.,Y,')) dr 

\ J s+a 

+ eE(ci>,(y)[(v,/(y/+j,G(x,v,,y/+j) 

-(v,/(y/+J,G(x,%,y/+j)]) 

+ e¥.Us{Y) r\vfiY,'),VyGiX',X)G{X',X))dr 

\ J s+a 

\ Js+a / 

We choose 5 > small enough, such that t + 5 <T, and deduce from the 
last identity that 



(26) 



k(<^s{Y) [f{Y,%J-f{Y^%J]da 

= e(^s{y) ['da r^"(v/(y/), 

V Jo Js+a 

y/) + VJ/G'(x^^ y/)G(x^^ y/)) dr 

+ Ef$,(y) ['da ['^\d^f{Y.nG{X^.X),G{X',,Y,'))dr 

\ Jo Js+a 

+ eE(^^s{Y) da[(V,/(y/+J,G(X,V„y/+J) 



-(v,/(y/+J,G(x,%,y/+j)] 
+ eEf$,(y) ['da r^"(v/(y/),Vj,G(x,%y/)G(x,%y/))dr 

V JO Js+a 

+ eEf^.,(y) ['da ^^^a2/(:^/)G(x,^y/),G(x,^y/))dry 

V JO Js+a J 

It follows from Lemma 5 in [11] that for any < s < t < T, 

Av/(y/), y/) + VJ,G(x,^ y/)G(x,^ y/) - 6(y/)) dr ^ 

J s 



and 



^Tra2/(>7)[G(X^^y/) ® G(X^^y/) - ia(y/)] dr ^ O 

J s 



in probability, as e — > 0. 

We can then take the limit in (26) as e — > 0, divide by 5 > 0, and let 5 — > 
since the process Y is right-continuous, yielding that for all / G G^(M^), 
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all < s < t < T and all bounded and 5-continuous functional defined 
on B([0,r];M^), 

E{[f{Yt)-f{Y,)mY)) 

(27) 

= e(^<!>s{Y) ^ [{Vf{Yr),b{Yr)) + \ Tr aV(>;Oa(>^r)] dr) , 
or in other words that 

Ml := /(Yi) - /(y,) - l\{Vf{Yr),biyr)) + i Trd^f{Yr)a{Yr)] dr 
is a martingale. 

It remains to show that t — > is a.s. continuous from [0, T] into M^, which 
is done in the following. 

Proposition 6. Let {y^, < t < T} he an I -dimensional semimartingale 
such that for aUl<i<i, all f £ C^iR), 

:= /(^/) - f{y^) - flfiYDUY,) + y"{Yl)a,,{Y,)] ds 
Jo 

is a martingale. Then {14,0 <t< T} is continuous. 

Proof. We note that the assumption implies that V / G C°°(M), Ml''^ is 
a local martingale. Hence in particular, for each 1 <i <i, 

Ml = Yl-Y^- fh{Y:)ds 
Jo 

is a local martingale, where Y^ denotes the ith component of Yt, and it 
follows from Ito's formula for (possibly discontinuous) semimartingales (see, 
e.g., [15], page 72) that V / G C°°(M), 

f{Y;) = f{Y^)+ f nYlMYs)ds+ f f{Y:-)dMl + \ ff"{Y:-)d[M% 
Jo Jo Jo 

+ E {f{Y:)-f{Y:.)-f'{Y:.)AY:-y'\Y:.){AY:f), 

0<s<t 

where {[M*]t,t > 0} denotes the quadratic variation process of the martin- 
gale M*. In the particular case f{y) = (y*)^, this identity reads 

(y/)2 = (yj)2 + 2 f' Y:k{Ys) ds + 2 f' y;_ dMi + [m%. 

Jo Jo 
Writing the assumption in the case f{y) = (j/*)^, we obtain that 

Mf := (y/)2 - (Y^)^ - [\2Y:k{Ys)+auiY,)]ds 

Jo 
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is a local martingale. Comparing the last two identities, we deduce that 
[M*]t — /q ajj(Y'<j) ds is a local martingale. Next, comparing the two different 
ways of writing (Y^)^ and using the identity 

iY:f = {y:_ f + 3(y;_ fAY: + syi {^Yif + {^Y^f, 

we deduce that Yl,o<s<t{^Y^)^ is a local martingale, from which we deduce, 
by comparing the two different ways of writing (1^*)^ ^^id using the identity 

{Y:f = {Yl_f + A{Y:_fAYl + Q{Y:_)\AYlf + Wi_{AYlf + (Ay;)^ 

that X]o<s<t('^^s*)^ ^ local martingale, which is impossible, unless it is 
identically zero. Since this result holds for any l<i<l, the proposition is 
established. □ 
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